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Abstract 

We present a small network for the testing of the entanglement of two 
ballistic electron waveguide qubits. The network produces different out- 
. put conditional on the presence or absence of entanglement. The structure 

' of the network allows for the determination of successful entanglement op- 

f^**) , erations through the measurement of the output of a single qubit. We also 

l/^ ■ present a simple model of a dynamic coulomb-like interaction and use it to 

' describe some characteristics of a proposed scheme for the entanglement 

of qubits in ballistic electron waveguides. 

£h ! 1 Introduction 

<ti , 

^ t A simple quantum computer consists of an array of qubits and a series of gates 

O^ 1 formed by single-qubit and two-qubit unitary transformations. A single qubit 

gate rotates the state of the qubit. The two-qubit gate creates an entangled 
pair of qubits. Any proposed system for quantum computation must provide a 
mechanism for pairwise entanglement of qubits. 

The possibility of performing quantum computation in ballistic electron 
waveguides was first proposed by Ioniciou et al 1 . In their approach, a sin- 
gle electron wavepacket and two parallel waveguides are used to form a "flying 
qubit," with one waveguide designated as the |0) state and the other as the |1) 
state. Subsequent work by Akguc et.al |2] and Snyder and Reichl 3 focused on 
the computation of stationary states of networks (rather than use time evolution 
of wavepackets) of such qubits. They showed that it is possible to obtain sta- 
tionary state solutions to the Schrodinger equation for fairly complex quantum 
networks of qubits and quqits. 

One proposed mechanism for entanglement of waveguide based qubits is the 
Coulomb interaction between electrons in different qubits. For example, a small 
segment of the waveguides in two qubits (which we call qubit A and qubit B) 
which represent the |1) state could be brought close to one another, or could 
be separated by a dielectric that allows interaction between electrons in the |1) 
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waveguides. This must be done in such a manner that electrons cannot tunnel 
between the qubits (see Fig. 1). If electrons in the |1) waveguides pass the 
interaction region at the same time they can interact and create a phase change 
in the network state |1,1) with no change in the remaining states |1,0), |0, 1) 
and |0, 0). This is sufficient to entangle the network. In Akguc et.al [2] a simple 
static model of this entanglement mechanism showed that a phase change of 
e m could be achieved for the state |1, 1). In subsequent sections, we analyse 
a dynamic model of the electron scattering process in the interaction region 
which confirms this prediction. We also analyse a simple two-qubit network 
which could allow a test for the efficiency of this entanglement mechanism. 

The waveguide structures we consider can be formed at the interface of a 
GaAs/AlGaAs semiconductor heterostructure. At temperatures, T ~ 0.1 — 
2.0 K, an electron travels ballistically with a phase coherence length of the 
order, ~ 30 — 40/im 0]. The gate structures themselves have been shown 
to be anywhere from 0.17^tm to 0.4^m in length Ej. The small network 
presented here contains few enough gates to be realizable with the coherence 
length presently achievable in semiconductor heterostructures. 

In Sect. II, we will present an electron waveguide network consisting of two 
qubits and a series of single-qubit and two-qubit gates which can be used to test 
for entanglement. We first construct the network with ideal single-qubit trans- 
formations and ideal entangling two-qubit transformations. We also construct 
the network with non-ideal entanglement gates and compare the outputs to the 
idealized case. We will see that the output of the non-ideal network could be 
used to determine whether or not the entanglement gates behave as expected. 
Then in Sect. Ill, we present a simple model of an electron scattering process 
that can achieve entanglement of a pair of qubits. A classical description of the 
dynamics is first discussed and then a steady state quantum scattering analysis 
of the same model is used to describe the behavior of the mutual phase acquired 
by the entangled electron current in the pair of waveguides. In Sect. IV, we 
make some concluding remarks. 

2 Entanglement Testing Network 

In this section, we describe a simple quantum network that can test the effec- 
tiveness of an entanglement gate. The network is shown in Fig. 2. It consists of 
a sequence of single qubit V NOT gates, Q, and two-qubit entanglement gates, 
V. (In Akguc et al it was shown that a single qubit V NOT could be con- 
structed in electron waveguides by using a properly constructed cavity which 
connects the two waveguide leads of the qubit.) Electrons are injected into the 
network from the left in a state = ci|l, 1) + c 2 |l,0) + C3|0, 1) + C4|0,0). 

This state is then acted on by a sequence of gates 

N = Q B -Qa.Vab.Q%.Q 3 a .Vab.Qb.Qa (1) 
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All three matrices act on the state vector $l = (ci, c 2 , C3, C4) T , where T 
denotes transpose. We write the two-qubit entanglement matrix in terms of 
phases fa , fa and so we can describe some general features of this matrix. 

When N acts on the input state we obtain an output state |<I>_r) = 

7V|$l) that gives the distribution of electrons exiting the quantum network on 
the right. For example, if <j)\ = <f>2 = and 6 — tt, an input state = 1 1 , 1) on 
the left leads to an output state |$r) = e JL 2 ZL |0, 1) on the right. In this example 
a series of single qubit operations and two entanglement operations produces 
an uncntangled output state. Although the output state is unentangled, the 
particular form of the output state will depend upon a successful entanglement 
of the qubits in the middle of the computation. 

We can indirectly test if the gate V is successful in entangling the two qubits 
by means of the network outlined above. A specific realization of the gate 
V is defined by the choice of the parameters <j>\, fa, and 9. Through these 
parameters we define two types of the gate V, one which entangles the qubits 
and one which does not. As the parameters are varied the output of the network 
is found for both types of V gate. We find that the entangling gate and the 
non-entangling gate give very different output in both the two qubit and one 
qubit bases, allowing for the determination of successful entangling operations 
in the network through the measurement of the output of only one of the qubits. 

A perfect entanglement gate V changes the phase only of the two-qubit 
state 1 1,1) and is represented by V where fa = fa = and 9 = n. A two- 
qubit gate that does not entangle the qubits changes the phase of the single 
qubit states and so that the two qubits remain separable. Such a 
gate is represented by the matrix V where fa + fa = 9. Due to the spatial 
symmetry of the quantum network we can choose fa = fa = 9/2 to represent 
a non-entangling two-qubit gate. Therefore, if we begin with an input state, 
= |1, 1) and act on it with a network, N, containing the perfect entangling 
gate V, where 9 = n and fa = fa = 0, wc find as above, N\l, 1) = e'^jO, 1). 
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A network containing the non-entangling gate where <f>\ = <f>2 = § and 9 = tt 
gives iVj 1, 1) = e l Tr|l,0), which is easily distinguishable from the case when 
entanglement is present. In Fig. 3, we plot the probability, P = |($#|0, 1)| 2 , of 
finding the ideal output \$>l) — |0, 1), as a function of 9 for both the entangled 
case and the unentangled case. We find that the respective outputs are most 
different when 9 = it, and equal when no phase change occurs. 

We can find the amount of probability exiting an individual waveguide in a 
given network by 

Prob A (\l)) = ($ R \(\l)(l\)\$ R ) (5) 

The probability of finding electrons in the |0) and |1) states of qubit A for both 
the entangled network and the unentangled network is plotted in Fig. 4. We 
see that the amount of probability exiting the |1) waveguide in relation to the 
|0) waveguide of qubit A is much greater for the entangled network than the 
unentangled network for phases angles near ir. There is a significant range of 
phase angles when the two networks would be distinguishable. 

In the situation above we have used the spatial symmetry of the network to 
set 0i — 4>2 —9/2 for a non-entangling gate V. An imperfect non-entangling 
gate need not split the phase angle 9 equally between the two individual qubits. 
We then write 4>2 = 9 — <j>\ and, given an input state = |1, 1), we analyze 
the output of the network as both 9 and (pi are varied. We find that when 
9 = 7r the output states of the entangling network and the non-entangling 
network are distinguishable for all values of 4>i. As above, the networks are 
most distinguishable when <p% = >p2 = tt/2. For the entangling network, the 
probability of finding the output state \$r) = |0, 1) when 6 = tt is 1. For the 
non-entangling network, the probability of finding the output state \$>r) — |0, 1) 
when 9 = n is never larger than 1 /4 for all values of <f>i . 

3 Dynamic Model of Coulomb Entangler 

In Akguc et.al 0, we presented a static model of Coulomb coupling between 
electrons in separate leads of a waveguide quantum network. We considered 
two parallel waveguide leads belonging to separate qubits, corresponding for 
example to the |1) states in the two qubits. We introduced a dielectric window 
between the leads that allowed electrons in the two leads to interact via their 
Coulomb interaction if they pass the dielectric window at the same time. We 
assumed that each electron produces a repulsive potential barrier in the path of 
the electron in the opposite lead. We then found that for certain energies the 
electrons they could resonantly pass the barrier and create a phase shift of e~ 
for each electron state giving an overall phase shift of e m for the network state 
|1, 1). In this and the next sections, we revisit that picture but with a dynamic 
model of the actual scattering process. 

In order to obtain an exactly soluble model of the scattering process, we 
simplify the model slightly. In the waveguide network, the actual scattering 
process takes place in the fixed (in space) dielectric window if two electrons (in 
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different waveguides) pass that window at the same time. The interaction they 
feel will be that of a finite range repulsive pulse (due to their mutual Coulomb 
interaction) whose shape, width and strength is determined the shape and width 
of the dielectric window and the distance between the waveguide leads. In our 
dynamic model we will neglect the dependence on the repulsive interaction due 
to the finite transverse width of the waveguide leads and we will allow the 
electrons to interact when they come within the range of their mutual repulsive 
interaction. We will choose our initial conditions so that this interaction occurs 
in a certain interval of space. Below we first consider a classical version of the 
model and then we consider the fully quantum scattering process. 



3.1 Classical Model of Coulomb Entangler 

Let us consider two one-dimensional straight wires, infinitely long in the x- 
direction, and separated by a distance, d in the y-direction. Electron A travels 
in the upper wire and electron B travels in the lower wire. Both electrons 
travel in the positive x-direction in their respective wires. We assume that 
the two electrons have nearly the same kinetic energy. Their velocities differ 
only by a small amount so that va = vq + Aw and vb = vq — Av. Initially 
the separation of the two particles in the x-direction is large enough that no 
appreciable interaction takes place. We assume that electron A is initially to 
the left of electron B but is closing the gap between them as they move up the 
x-axis. 

We can write the total Hamiltonian for the system in the form 

H = -^-Va + ^-Pb 2 + TJ-. ^ 77=E tot , (6) 

2m 2m cosh z (a(x A - x B )) 

where pa — vuva and xa (pb — tuvb and xb) are the momentum and position 
of particle A (particle B), Vb is the maximum interaction strength, — is the 
width of the interaction potential between the two electrons and E tot is the 
total energy of the system. In a Coulomb-like interaction the distance between 
the wires, d, determines the maximum interaction strength, Vb, between the 
particles but otherwise does not add to an understanding of the interaction 
itself. 

The center of mass momentum and position of the electrons are P = pa+Pb 
and X = + Xb), respectively. Their relative momentum and position are 

x = xa — x B and p = ^(pa — Pb), respectively. In terms of these coordinates, 
the Hamiltonian takes the form 



4m m cosh 2 (era) 



H =— + -+ J, - =Etof (7) 



We see that the center of mass momentum and the center of mass energy E cm = 
P 2 1 Am are constants of the motion. 

All the interesting dynamics occurs in the relative motion of the two electrons 
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whose Hamiltonian is given by 



Vo 

cosh 2 (aiE) 



E r 



(8) 



where E r is the energy contained in the relative motion of the particles. The 
character of this motion is determined by the relationship between the energy 
of relative motion, E r = E to t — E cm , and the interaction strength, Vq. The 
phase space diagram for the relative motion is plotted in Fig 5. Electrons with 
relative energy < E r < Vo interchange their momenta during the collision but 
not their relative postions (the phase space motion corresponds to the curves 
that cross the x-axis in Fig. 5). Trajectories with relative energy Vo < E r < oo 
interchange their position and not their relative momenta during the collision 
(the curves that cross the p-axis in Fig. 5). 

The case where electrons A and B have approximately the same velocity so 
va = v o + Av and vb = vq — Av with Ai><Cvo, and both travel in the positive 
x-direction corresponds to the case < E r < Vq. Electron A will catch up to 
electron B and they will undergo a collision with the result that they interchange 
their velocities but not their positions. We combine the solutions for the center 
of mass coordinate and the relative coordinate for the case E r < Vq and obtain 



x A (t) 



m 2a 



— — 1 cosh [ — a\ — -t 
E r \ V m , 



(9) 



and 



x B (t) = 



m 2a 




(10) 



For these solutions, the interaction is centered at x = at time t=0. In the 
asymptotic regions where both particles are far away from the interaction (t — > 
oo, t — > — oo) the particles move with constant velocity. The particles exchange 
velocity during the interaction and do not pass each other. 

In ballistic electron waveguides built using GaAs-AlGaAs heterostructures 
the energy of the traveling electrons at low temperatures is very close to the 
Fermi energy of the electron gas |B] . We would therefore expect that the energies 
of any two electrons traveling through a coulomb coupler-like structure would be 
quite similar, resulting in a small relative energy with respect to the interaction 
potential. This corresponds classically to the case E r < Vq considered above. 



3.2 Quantum Scattering Model of Coulomb Entangler 

Let us now consider the quantum realization of the classical model described 
above. The Schrodinger equation for the two particle system is 



where ^ = ^(xa,xb) is the energy eigenstate of the two particle system. If 
we again change to center of mass and relative coordinates, the Schrodinger 
equation takes the form 

n 2 (i d 2 d 2 \ v y 

where VP is now a function of the center of mass and relative coordinates. 
The center of mass momentum and energy are again constants of motion for 
this system. We assume a separable form for the two-particle wavefunction, 
ty(X,x) — tp(X)(j)(x). The solution for the the center of mass wave function is 

4>(X) = e lKX (13) 

where K = P/% — J AmE cm /h 2 is the center of mass wavevector. The solution 
for the wavefunction describing the relative motion is 



-ik —ik —ik 1 , 

8, +1,-1-0 

a a a 2 



(14) 



| -1 + Jl - AmV /a 2 h 2 . The center of mass solution, i/j(X), is chosen to 



where F is a hypergeometric function, k = yj mE r /h 2 , ( — tanh(ax) 



2 

represent the state of two particles traveling in the positive x-direction and is 
normalized to unity. 

We are considering the scattering of two electrons, for the case E r < Vq, 
traveling along the pair of waveguides in the positive x-direction. Initially the 
electrons enter the system from the left such that electron A begins to the left 
of electron B and electron A has a slightly larger momentum than electron B. 
The repulsive interaction potential between the two electrons falls off rapidly 
enough that asymptotically (£— >±oo) the electrons are free. 

From the discussion of the classical version of this problem (for E r < V ) we 
see that there are two asymptotic regimes. In one regime (x— > — oo), electron 
A remains to the left of electron B, but during the collision they interchange 
momenta (this is the only case that is allowed classically). However, quantum 
mechanically the regime (x— > + oo) is also allowed. This would require the wave- 
function to tunnel through the barrier in the relative motion problem. We can 
now write the asymptotic form of the solution for the relative motion problem 
in the form 

<p(x -»■ oo) = Te tkx and <f>(x -» -oo) = e lkx + Re~ ikx . (15) 

The coefficient T is the probability amplitude that the electrons interchange 
position and not momentum during the collision. The coefficient R is the prob- 
ability amplitude that the electrons interchange momentum and not postion 
during the collision (the classically allowed case). The term e lkx is the wave- 
function for the relative motion before the collision. If we take the asymptotic 
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limits (x— >±oo) of the hypergeometric function we obtain the following expres- 
sions for the probability amplitudes T and R 



r(=ik- s )T(=i!z+s + l) 



T = , — Mr s (16) 



r( f)r(f -,)r(f + s + i) 
r(^)r(- s )r( s + i) 



where T(x) is the gamma function. 

We can now write the total wavefunction for the system in the asymptotic 
regions {xa~ > — oo,xs^ — oo) and (xa^ + oo,xb^ + oo). For (x^^ — oo,xs^ — 
oo ) the total wavefunction is 

^{x A ,x B ) =e ikAXA e ikBXB , (18) 

where ItA = pa/^ and ks = Pb/Ti are the incident wavevectors of electrons A 
and B. For (x^— > + oo,xb — > + oo) the total wavefunction is 

*(x A , x B ) = Te ikAXA e ikBXB + R e lk B*A e ik A x B ^ 

This simple model of Coulomb entanglement predicts no reflection of an 
individual electron due to the collision, but simply a mutual phase shift of the 
two electrons and a possible exchange of momenta. This bodes well for future 
implementations of such structures in quantum processing devices, as reflection 
of individual electron probability at the computational gates plays a large role 
in determining the fidelity of a computation 13] . 

As stated above, in ballistic electron waveguides in GaAs-AlGaAs semi- 
conductor heterostructures the incoming energies of each electron is expected 
to be near the Fermi energy of the device, Ef. We assume that widths of 
the waveguides are equal so that the energy required for the first transverse 
mode is the same. For electrons in the first propagating channel of the waveg- 
uide leads, this means that the momenta of the electrons will be given by 

k A = y ^W) _ an d k B = ^"(W) I (If where S E is deviation 
in energy from the Fermi energy due to the finite temperature of the semicon- 
ductor material. For low temperatures we can expect 8e to by very small and 
therefore the relative momentum of the two electrons to be very small. From 
above discussion, we see that as the relative momentum k — > 0, R — > — 1 and 
T — > 0. This would correspond to the electrons exchanging momentum and 
leaving the interaction region with a mutual phase change of e 47r . This is just 
what is needed to obtain optimum entanglement in the network described in 
the previous section. 

To determine how the reflection and transmission amplitudes might vary in 
an implementation of the network at finite temperature, we use numerical values 
similar to those use in GaAs-AlGaAs quantum networks. We define a unit 
of length, luq — AOnm and a unit of energy, 

£ = t-^-tt = 0.000355ey, where 
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m = 0.067m e is the effective eiectron mass in GaAs-AlGaAs semiconductor 
structures and m e is the mass of the free electron. If we assume that the leads 
have a transverse width w — 160A, the the electrons propagate in the first 
channel for Fermi energies 61.7<Ef / Eq<246.8. We use an interaction potential 
of Vq — 32.14Ex) and an interaction region of length 1/a = Fig 6 shows the 
behavior of the phase angle of the reflection amplitude as the relative momentum 
varies. For small values of the relative momentum the reflection amplitude, 
and therefore the mutual phase between the electron currents, approaches e m . 
Fig 7 shows how the reflection probability varies with the relative momentum. 
Reflection dominates for small relative momentum. 

In order to maintain quantum coherence in these types of devices, temper- 
atures must be on the order of a few degrees KelvinJHj. When both qubits 
are formed at the same semiconductor heterostructure we can assume the same 
Fermi energy value in both qubit structures. At the low temperatures associated 
with these types of semiconductor devices the electrons travel with an energy 
very near the Fermi level. Therefore the average deviation from the Fermi level 
of each of our two traveling electrons corresponds to the relative energy of the 
electrons incident on the interaction region. At small temperatures the electron 
energy deviates from the Fermi energy an average amount Se ~ ksT, where 
ks is Boltzmann's constant and T is the temperature. For a temperature of 
AK we find the average separation in energy to be Se ~ Eq. This allows us to 
find the average deviation in each electrons longitudinal momentum, 5k- Taking 
values for the leads of width w = 160A and Fermi energy Ef /Eq — 150 we find 
5k ~ 0.0027rtm~ 1 . Then using a relative momentum of k — 5k we find near unit 
probability for momentum exchange and an amplitude phase angle very near tt 
so that the phase angle of the amplitude, R = e 10 , is 9 = tt± 5$ where 5$ = 0.13 
radians. 

4 Conclusions 

We have presented a network for the testing of entanglement in ballistic elec- 
tron waveguide qubits. The entangling properties of the coulomb gate are dis- 
tinguishable for phase angles close to tt. The simple model of a coulomb-like 
coupler predicts a mutual phase angle of the |1, 1) state very near tt when the 
relative momentum between the two particles is very small. There is no re- 
flection of individual electrons at the coulomb region. All incoming probability 
continues forward through the coulomb coupler region towards the output side 
of the network. 
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Figure 1: Each qubit is a pair of waveguides. The spatial location of the electron 
in the waveguides determines the state of the qubit. Here both qubits are in 
state |1). The waveguides representing state |1) are brought near each other to 
facilitate the coulomb interaction of the electrons, effecting a two-qubit unitary 
transformation. 



Figure 2: A schematic of the entanglement testing network. Boxes represent 
individual transformations. Q and Q 3 are single-qubit transformations. V is a 
two-qubit transformation. 



Figure 3: A plot of the probability of finding the output state |1, 0) as a function 
of the phase angle 9 for both the entangled and unentangled cases, (a) Entangled 
network with </>i =02 =0. (b) Unentangled network with <pi = <p 2 = 9/2. 
When 9 is near 7r the entangled and unentangled situations are most easily 
distinguished. 
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Figure 4: The output probability for each state of qubit A. (a) Entangled net- 
work with fa =02 = 0. (b) Unentangled network with <j>i = <p2 = 9/2. For 
the entangled network the amount of electron current in the \1)a state is large 
compared to the |0)a state for phase angle 9 near n. The opposite is true for 
the unentangled network. 



Figure 5: The phase space of the relative motion. The dashed line is the sep- 
aratrix between electrons with relative energy < E r < Vo and electrons with 
relative energy Vo < E r < oo. Electrons with relative energy < E r < Vo 
(inside the separatrix) interchange their momenta during the collision and not 
their relative postions. Trajectories with relative energy V < E r < oo (outside 
the separatrix) interchange their position and not their momenta during the 
collision. (All units dimensionlcss.) 



Figure 6: Plot of the phase angle 9 of the probability amplitude R — e 10 ver- 
sus relative momentum, k. The phase angle is very near n when the relative 
momentum is very small. 



Figure 7: The probability |i?| 2 of momentum exchange during the collision. 
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